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Abstract 

The  circulation  sensitivity  to  the  choice  of  wind-forcing  product  is  investigated  with  the  NRL  Layered  Ocean  Model 
(NLOM)  configured  for  the  Japan/East  Sea.  Monthly  climatologies  from  seven  different  wind-stress  data  sets  (and 
wind-stress  curl)  are  formed  from  observed  and  model-derived  atmospheric  data  sets.  The  resulting  maps  of  wind-stress 
curl  reveal  significantly  different  spatial  patterns  and  magnitudes,  even  two  with  largely  opposite  large-scale  patterns  of 
wind-stress  curl.  These  wind  sets  were  used  in  forcing  three  sets  of  simulations,  1/8°  linear  and  1/8°  and  1/32°  nonlinear. 
In  addition,  seasonally  varying  straits  forcing  (inflow  through  Tsushima  balanced  by  outflow  through  Tsugaru  and 
Soya)  was  included  in  all  the  simulations,  and  simulations  with  straits  forcing  only  were  performed.  The  1. 5-layer  linear 
reduced-gravity  simulations  include  only  the  lowest  order  dynamics,  mainly  Munk  (3 1/3  western  boundary  layers  (due  to 
both  wind  and  straits  forcing)  and  a  Sverdrup  interior.  The  nonlinear  simulations  add  bottom  topography,  multiple 
internal  modes,  diapycnal  mixing,  and  ventilation  of  layer  interfaces.  At  1/8°  resolution,  only  weak  barotropic/ 
baroclinic  instabilities  occur,  but  at  1/32°  resolution  these  are  much  stronger,  allowing  vigorous  transfer  of  energy  from 
the  upper  ocean  to  the  abyssal  layer  via  baroclinic  instability.  This  drives  much  stronger  mean  abyssal  currents  that 
more  strongly  steer  upper-ocean  current  pathways  than  at  1/8°,  i.e.  there  is  much  stronger  upper  ocean-topographical 
coupling.  The  linear  model  simulates  most  of  the  basic  features,  e.g.,  the  subpolar  gyre  with  all  but  the  straits  forcing 
only,  the  East  Korean  Warm  Current  (EKWC)  and  its  connection  to  the  subpolar  front  with  all  but  one  wind-forcing 
set,  but  the  remaining  wind  set  gives  a  continuous  Nearshore  Branch  (NB)  of  the  Tsushima  Warm  Current  along  the 
coast  of  Honshu.  In  all  of  the  linear  simulations  with  an  EKWC,  the  separation  latitude  from  the  coast  of  Korea  is 
determined  by  the  latitude  (~41°N)  of  the  Tsugaru  Strait  outflow  (too  far  north)  except  for  one,  which  is  also  too  far 
north.  All  of  the  nonlinear  simulations  include  the  NB,  which  is  mainly  due  to  isopycnal  outcropping.  This  creates 
another  inflow  to  outflow  pathway,  diverting  some  of  the  flow  from  the  EKWC  and  allowing  the  wind-stress  curl  to 
compete  better  with  straits  forcing  in  determining  the  EKWC  separation  latitude.  Using  the  1/8°  nonlinear  model,  three 
of  the  wind-forcing  sets  yield  a  realistic  separation  latitude  for  the  EKWC  (~37-38°N)  and  two  more  are  marginally 
successful.  The  two  remaining  and  the  simulation  with  straits  forcing  only  yield  EKWC  overshoot  to  ~40°N.  At  1/32° 
resolution,  one  of  the  simulations  (one  with  forcing  that  gave  marginally  realistic  separation  at  1/8°)  has  a  realistic 
EKWC  separation  latitude.  The  three  wind  products  that  yielded  realistic  EKWC  separation  at  1/8°  continue  to  give 
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realistic  separation  at  1/32°  via  the  wind-stress  curl,  but  the  remainder  give  realistic  separation  via  upper 
ocean-topographical  coupling.  At  1/32°  resolution  all  but  one  of  the  simulations  showed  a  strong  Yamato  Rise 
impact  on  the  pathway  of  the  subpolar  front  due  to  steering  by  abyssal  currents.  At  1/8°  resolution  only  two  of  the  eight 
showed  any  such  impact,  one  due  to  abyssal  currents  driven  by  vertical  mixing  in  the  subpolar  gyre.  In  all  the  other  1/8° 
simulations  the  mean  abyssal  currents  were  weak,  but  in  one  case  weak  abyssal  currents  were  aided  by  a  region  of 
strong  negative  wind-stress  curl  not  seen  in  the  other  wind  products.  In  addition,  regions  of  the  ventilation  of 
intermediate  and  deep  water  in  response  to  each  wind  set  are  identified.  A  series  of  six  1/32°  simulations  forced 
interannually  with  6-hourly  winds  1979-93  are  used  to  separate  the  model  variability  due  to  wind  forcing  (deterministic) 
and  flow  instabilities  (nondeterministic).  This  analysis  indicates  that  flow  instabilities  are  the  dominant  source  of  SSH 
variability  over  a  large  majority  of  the  basin,  especially  near  37-38°N,  where  the  EKWC  separates  from  the  coast,  and 
along  the  pathway  of  the  subpolar  front  following  the  north  rim  of  the  Yamato  Rise  towards  the  Tsugaru  Strait.  Wind 
forcing  dominates  the  SSH  variability  in  the  coastal  wave  guide  and  over  the  deep  basins  and  troughs. 

Published  by  Elsevier  Ltd. 


1.  Introduction 

The  circulation  in  the  Japan/East  Sea  (JES, 
hereafter)  can  be  significantly  influenced  by  many 
physical  factors,  including  inflow  through  the 
Tsushima/Korea  Strait  (with  accompanying  out¬ 
flow  through  Tsugaru  and  Soya  Straits),  surface 
wind  forcing,  surface  buoyancy  forcing  (heating 
and  cooling,  salt-  and  fresh-water  fluxes),  and  the 
bottom  topogaphy.  There  are  also  dynamical 
factors  that  affect  the  circulation,  such  as  iso- 
pycnal  outcropping  and  barotropic  and  baroclinic 
flow  instabilities.  These  factors  impact  the  surface 
and  deep  circulation  on  different  space  and  time 
scales,  ranging  from  sub-mesoscale  to  interannual. 
All  are  believed  to  have  a  significant  impact  on 
the  location  and  strength  of  the  large-scale 
circulation  features.  In  recent  years  much  greater 
knowledge  and  understanding  of  the  mean  current 
structures  in  the  JES  have  evolved,  but  so  has  the 
recognition  that  there  is  a  rich  mesoscale  eddy  field 
and  that  even  the  large-scale  current  systems 
exhibit  significant  variability  over  a  wide  range 
of  frequencies.  However,  several  schematics  of 
the  general  JES  circulation  provide  a  guide  for 
evaluating  JES  modeling  efforts.  The  schematic 
shown  in  Fig.  1  depicts  the  mean  surface  circula¬ 
tion  features,  including  the  East  Korean  Warm 
Current  (EKWC)  (a  branch  of  the  Tsushima 
Warm  Current  (TWC)),  the  Liman  and  North 
Korea  Cold  Currents  (LCC  and  NKCC,  respec¬ 
tively),  the  subpolar  front,  and  the  Nearshore 
Branch  (NB)  of  the  TWC.  Other  schematics  of  the 


JES  show  slight  variations  or  more  detail,  as 
shown  by  Naganuma  (1977)  and  Yarachin  (1980), 
but  the  large-scale  features  are  similar. 


130°E 


TSS  :  Tsushima  Strait 
TGS  :  Tsugaru  Strait 
SS  :  Soya  Strait 
MS  :  Mamiya  Strait 
YR  :  Yamato  Rise 
JB  :  Japan  Basin 
YB  :  Yamato  Basin 
UB  :  Ulleung  Basin 
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Fig.  1  Schematic  map  for  sea-surface  currents  in  the  JES. 
Abbreviations  explained  within  figure.  From  Morimoto  and 
Yanagi  (2001),  after  Senjyu  (1999). 
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To  lowest  order,  ocean  dynamics  are  largely 
governed  by  the  wind-stress  curl  through  the 
Sverdrup  relation.  Ocean  circulation  models,  there¬ 
fore,  require  the  surface  wind-stress  or  its  curl  as  a 
specified  boundary  condition,  and  climatological 
fields  are  commonly  used.  Several  modeling  studies 
in  the  JES  have  been  undertaken  in  the  last  two 
decades,  including  pioneering  studies  by  Yoon 
(1982)  and  Kawabe  (1982),  and  more  recently 
studies  by  Sekine  (1991),  Seung  and  Kim  (1993), 
Holloway  et  al.  (1995),  Kim  and  Yoon  (1999),  and 
Hogan  and  Hurlburt  (2000).  These  modeling 
studies  yield  a  wide  range  of  numerical  solutions 
that  can  be  compared  to  observations  or  schematic 
diagrams  such  as  Fig.  1  to  assess  their  perfor¬ 
mance,  but  the  variety  of  ocean  model  designs, 
configurations,  and  forcing  functions  can  make  it 
difficult  to  distill  the  source  of  the  differences  and 
to  determine  which  have  the  greatest  impact. 

In  this  paper,  all  of  the  simulations  include  the 
throughflow  forcing  via  the  straits.  Given  this,  we 
show  that  the  choice  among  seven  different 
monthly  wind-forcing  products  can  have  a  pro¬ 
found  impact  on  numerical  ocean  model  simula¬ 
tions  of  all  the  current  systems  shown  in  Fig.  1  plus 
many  other  features  of  the  circulation,  including 
the  abyssal  circulation.  For  example,  some  give 
realistic  separation  of  the  EKWC  from  the  east 
coast  of  Korea,  while  others  do  not.  In  addition  to 
wind-stress  curl,  Hogan  and  Hurlburt  (2000) 
(hereafter  HH2000)  discovered  another  mechan¬ 
ism  for  realistic  separation  of  the  EKWC,  upper 
ocean-topographical  coupling  via  baroclinic  in¬ 
stability.  Either  mechanism  by  itself  can  be 
sufficient  to  give  realistic  EKWC  separation  and 
to  impact  many  other  features  of  the  circulation, 
but  HH2000  needed  1/32°  resolution  for  each 
model  variable  to  obtain  EKWC  separation  via 
the  upper  ocean-topographical  coupling  mechan¬ 
ism  alone  and  for  this  mechanism  to  have  a  major 
impact  on  the  overall  circulation. 

In  this  paper,  the  role  of  wind  forcing  in  the 
absence  and  presence  of  upper  ocean-topographi¬ 
cal  coupling  is  investigated  by  performing  simula¬ 
tions  with  the  Naval  Research  Laboratory  (NRL) 
Layered  Ocean  Model  (NLOM)  at  1/8°  and  1/32° 
resolutions,  respectively.  Specifically,  at  1/8°  re¬ 
solution  we  have  performed  two  sets  of  eight 


simulations  (seven  simulations  for  each  of  the 
monthly  wind-stress  climatologies  and  one  with  no 
wind  forcing).  The  first  set  is  linear,  1.5-layer 
reduced  gravity,  and  includes  only  the  lowest  order 
dynamics,  basically  a  Sverdrup  (1947)  interior  with 
Munk  (1950)  (31/3  western  boundary  layers  due  to 
both  wind  and  throughflow  forcing.  The  second 
set  is  moderately  nonlinear  and  adds  bottom 
topography,  multiple  internal  modes,  diapycnal 
mixing  and  ventilation  of  the  Lagrangian  model 
layers  (buoyancy  forcing  is  neglected).  These 
additions  allow  the  possibility  of  barotropic/ 
baroclinic  flow  instabilities,  the  formation  of 
inertial  jets  and  boundary  layers,  and  isopycnal 
outcropping.  A  similar  set  of  simulations  was 
performed  at  1/32°,  except  that  the  simulations 
forced  by  the  Comprehensive  Ocean-Atmosphere 
Data  Set  (COADS)  and  throughflow  only  are 
omitted  due  to  the  similarity  of  the  simulations 
with  the  HR-forced  simulation  at  1/8°  resolution. 
The  impact  of  increased  horizontal  grid  resolution 
is  discussed  in  Section  6. 


2.  The  NRL  Layered  Ocean  Model 

NLOM  is  a  semi-implicit  primitive-equation 
ocean  model  that  casts  the  layered  model  equa¬ 
tions  in  transport  form.  The  model  retains  the  free 
surface  and  can  include  realistic  bottom  topogra¬ 
phy  as  long  as  it  is  restricted  to  the  lowest  layer. 
The  original  model  design  is  described  in  Hurlburt 
and  Thompson  (1980),  and  significant  enhance¬ 
ments  are  described  in  Wallcraft  (1991),  Wallcraft 
and  Moore  (1997),  and  Moore  and  Wallcraft 
(1998),  including  an  implementation  in  spherical 
coordinates. 

All  of  the  nonlinear  simulations  used  in  this 
study  have  four  Lagrangian  layers  in  the  vertical 
with  interfaces  initially  situated  at  60,  130,  and 
250  m  with  density  differences  (<7©)  across  each 
interface  of  0.80,  0.32,  and  0.42,  respectively. 
The  comparisons  using  different  atmospheric 
forcings  are  done  at  both  1  /8°  and  1  /32°  horizontal 
resolutions  between  like  variables  on  a  C-grid.  The 
role  of  the  1/8°  simulations  is  to  limit  the  impact 
of  flow  instabilities  while  maintaining  sufficient 
nonlinearity  to  assess  its  impact  on  the  response  of 
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current  systems  to  atmospheric  forcing.  The  1/32° 
simulations  include  the  impact  of  the  different 
wind  forcing  as  well  as  higher-order  dynamics 
(such  as  upper  ocean-topographical  coupling),  as 
described  in  HH2000.  An  additional  suite  of 
simulations  with  1/32°  resolution  is  used  to  assess 
the  impact  of  flow  instabilities,  and  wind/straits 
forcing  on  the  variability  is  described  in  Section  7. 
Aside  from  the  grid  resolution  (1/8°  vs.  1/32°),  the 
only  other  physical  (vs.  dynamical)  difference  is 
the  horizontal  eddy  viscosity  (50  vs.  5m2/s). 
Although  thermodynamic  versions  of  the  model 
exist,  all  the  simulations  used  here  are  purely 
hydrodynamic,  i.e.  they  contain  no  horizontally 
varying  density  gradients  within  a  layer.  Holloway 
et  al.  (1995)  found  buoyancy  forcing  to  be  less 
significant  relative  to  wind  forcing,  at  least  for  the 
large-scale  surface  circulation.  That  claim  is  not 
made  here,  since  all  of  the  simulations  in  this  study 
exclude  buoyancy  forcing;  rather  the  focus  of  this 
paper  is  the  response  to  the  wind  forcing. 

For  all  simulations,  the  throughfiow  is  provided 
via  open  boundary  conditions  at  the  Tsushima/ 
Korea,  Tsugaru,  and  Soya  Straits.  Inflow’  through 
the  Tsushima/Korea  Strait  is  divided  into  two 
separate  inflow  ports,  one  each  for  the  western  and 
eastern  channels  of  the  strait.  Consistent  with 
observations  (see  Preller  and  Hogan,  1998),  the 
mean  volume  transport  of  the  inflow  is  2.0  Sv 
(Sv  =  106  m3/s),  but  a  seasonal  signal  is  super¬ 
imposed  such  that  the  maximum  inflow  and 
outflow  of  2.66  Sv  occur  in  July  and  the  minimum 
of  1.34Sv  occurs  in  January  (33%  peak  deviation 
from  the  mean).  The  inflow  is  consistently  segmen¬ 
ted  so  that  75%  enters  through  the  western  channel 
and  25%  through  the  eastern  channel.  Vertically, 
two-thirds  of  the  inflow  (for  both  channels)  enter 
the  basin  through  layer  1,  while  the  remaining  one- 
third  enters  through  layer  2.  The  straits  are  closed 
in  layers  3  and  4.  Inflow  through  the  Tsushima/ 
Korea  Strait  is  exactly  balanced  by  outflow 
through  the  Tsugaru  and  Soya  straits  using  a 
modified  Orlanski  (1976)  boundary  condition.  The 
outflow  is  distributed  such  that  two-thirds  of  the 
flow  exit  through  Tsugaru  Strait  and  the  remaining 
one-third  through  Soya  Strait.  The  vertical  dis¬ 
tribution  for  the  outflow  is  the  same  as  for  the 
inflow.  HH2000  contains  a  detailed  description  of 


the  circulation  dynamics,  including  the  coupling 
of  the  abyssal  and  surface  circulation  and  the 
impact  of  bottom  topography.  It  also  includes  a 
complete  description  of  the  JES  model  configura¬ 
tion  used  here. 


3.  The  atmospheric  forcing  products 

The  surface  wind  stress  is  the  primary  source  of 
momentum  flux  into  the  ocean  and  the  primary 
forcing  of  the  upper-ocean  circulation.  Ocean 
circulation  models,  therefore,  rely  on  the  surface 
wind-stress  (primitive  equation  models)  or  its  curl 
(quasi-geostrophic  models)  as  a  specified  surface 
boundary  condition.  In  particular,  the  curl  of  the 
wind  stress,  which  is  important  dynamically 
because  it  is  a  source  of  vorticity,  is  the  principal 
forcing  mechanism  for  the  wind-driven  circulation. 
In  NLOM,  the  wind  forcing  is  input  via  wind- 
stress  and  prescribed  as  a  body  force  acting  only 
on  the  layer  in  direct  contact  with  the  atmosphere, 
although  the  top  layer  responds  to  the  curl 
through  Ekman  dynamics  and  the  Sverdrup 
relation. 

There  are  several  different  sources  of  wind-stress 
data  sets,  ranging  from  climatologies  formed  via 
statistical  analyses  of  marine  wind  observations  to 
those  formed  by  data-assimilative  operational 
atmospheric  models.  Hence  these  wind-stress  data 
sets  can  be  significantly  different  from  each  other 
and  subject  to  several  possible  sources  of  error. 
Those  constructed  from  historical  records  of 
marine  observations  are  particularly  vulnerable 
to  errors  based  on  the  estimation  of  wind  speed 
derived  from  ship  observations  at  sea  and  the 
parameterization  of  the  surface  momentum  flux  by 
a  bulk  aerodynamic  formula.  Model-based  data 
sets  also  are  influenced  by  grid  resolution,  model 
physics,  data  assimilation  methods,  and  first-guess 
errors.  One  shortcoming  both  wind-forcing  data 
sources  share  is  the  native  horizontal  grid  resolu¬ 
tion,  which  is  typically  much  larger  than  the 
internal  radius  of  deformation  for  the  ocean.  This 
is  particularly  exacerbated  in  the  JES,  where  some 
data  sets  with  relatively  coarse  resolution  (2.5°  for 
example)  may  have  only  a  few  grid  points  over  the 
JES.  As  expected,  OGCM  response  can  be  quite 
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big.  2.  Local  ion  map  of  open  boundaries  and  annual  mean  wind-stress  curl  in  10s  Pa/M  for  the  seven  different  wind- forcing  data  sets 
used  in  the  study,  Labels  indicate  which  product  is  depicted. 


different  depending  on  which  wind-forcing  data 
set  is  employed.  The  object  of  this  paper  is  to 
examine  and  quantify  the  differences  in  JES 
simulations  that  arise  from  forcing  an  ocean 
model  with  seven  different  wind-forcing  data  sets, 
which  arc  described  below  and  depicted  in  Fig.  2. 

The  Hcllcrman  and  Rosenstein  (1983)  (HR, 
hereafter)  monthly  mean  wind-stress  climatology 
has  been  widely  used  to  force  global  and  regional 
OGCMs.  This  climatology  is  comprised  of  106 
years  of  surface  marine  data,  from  1870  through 
1976,  derived  using  a  bulk  aerodynamic  formula, 
the  World  Meterological  Organization  (WMO) 
Beaufort  equivalent  scale  for  wind  speed,  and 
paramctcrizations  similar  to  Bunker  (1976).  The 
native  grid  resolution  is  2  x2c.  The  wind  stress 
was  smoothed  with  one  pass  of  a  nine-point 
smoother  to  eliminate  a  two-grid  interval  noise 
that  otherwise  gives  a  noodling  pattern  in  the  wind 


stress  curl  when  the  HR  wind  stress  is  interpolated 
to  the  1/8C  model  grid.  After  the  smoothing  the 
annual  mean  wind-stress  curl  (Fig.  2b)  shows  little 
small-scale  structure,  and  is  characterized  by 
positive  (negative)  wind-stress  curl  in  the  north¬ 
eastern  (southwestern)  part  of  the  basin. 

da  Silva  et  al.  (1994a,  b)  developed  a  new 
Beaufort  equivalent  scale  that  they  utilized  to 
improve  the  estimated  wind  speeds  of  each  observa¬ 
tion  (mostly  wind  estimates  over  1945-89)  in  the 
COADS.  Hence,  this  data  set  is  different  from  that 
originally  described  by  Woodruf  et  al.  (1987).  The 
new  analysis  also  used  higher  resolution  (F  x  I 
globally  excluding  Antarctica)  and  a  stability- 
dependent  value  for  the  drag  coefficient  (CD). 

Na  et  al.  (1992,  1999)  calculated  sea-surface 
winds  (at  10  m  above  sea  level)  twice  daily  for  an 
18-year  period  (1978-95)  utilizing  the  Cardone 
(1969)  model  and  twice-daily  weather  maps 
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published  by  the  Japan  Meteorological  Agency.  A 
bulk  formula  was  used  to  calculate  the  wind-stress 
using  the  drag  coefficient  described  by  Garratt 
(1977).  From  this,  monthly  mean  wind-stress  and 
wind-stress  curl  were  calculated  over  the  East 
Asian  Seas  region  with  0.5  x0.5c  resolution.  In 
the  JES  region,  the  annual  mean  wind-stress  curl 
(Fig.  2d)  is  characterized  by  positive  wind-stress 
curl  in  the  northern  part  of  the  basin  with  local 
positive  maxima  near  Vladivostok  and  the  East 
Korea  Bay,  while  negative  curl  is  present  in  the 
southeastern  part  of  the  basin  along  the  coast  of 
Honshu. 

The  European  Centre  for  Medium-Range 
Weather  Forecasts  (ECMWF)  routinely  performs 
global  operational  atmospheric  nowcasts  using  a 
data-assimilative  atmospheric  model.  An  ECMWF 
re-analysis  of  archived  data  was  used  to  form  two 
monthly  wind-stress  climatologies,  one  based  on 
the  10 m  winds  (EC10M,  hereafter)  and  the  other 
based  on  the  1 000  mb  winds  (EC 1 000,  hereafter). 
The  ECMWF  reanalysis  was  performed  at  6- 
hourly  intervals  over  1979  93,  using  a  fixed  1994 
version  of  the  ECMWF  operational  forecast  model 
and  data  assimilation  scheme  (Gibson  et  al.,  1997). 
The  reanalyses  have  the  advantage  over  the 
archived  operational  winds  in  that  observations 
not  available  in  the  original  nowcasts  are  included, 
and  the  entire  time  series  is  free  of  spurious 
variations  associated  with  real-time  changes  in 
the  operational  model  and  data  assimilation 
scheme  (although  changes  in  the  data  stream 
remain).  The  resolution  of  the  EC10M  winds  was 
spectral  T106/N80  Gaussian  output  on  a  1.125 
grid  (lat.  x  long.),  but  the  EC  1000  winds  were 
output  on  a  2.5C  x  2.5'  grid  after  truncation  to  T47. 
There  are  significant  differences  in  the  resulting 
wind-stress  curl  fields  (Figs.  2e  f),  which  are  due  in 
part  to  the  difference  in  horizontal  grid  resolution. 
The  EC10M  curl  shows  much  finer  scale  structure 
depicting  alternating  areas  of  positive  and  negative 
curl,  whereas  the  EC1000  curl  basically  shows 
positive  (negative)  values  in  the  northeastern 
(southwestern)  part  of  the  basin. 

A  monthly  mean  climatology  of  surface  stresses 
covering  July  1990-- June  1997  was  created  using 
archived  global  operational  atmospheric  model 
nowcasts  from  the  Fleet  Numerical  Meteorology 


and  Oceanography  Center's  Navy  Operational 
Global  Atmospheric  Prediction  System  (NO¬ 
GAPS)  (Hogan  and  Rosmond,  1991).  The  grid 
resolution  of  these  analyses  is  1.25'  x  1.25  .  This 
wind  forcing  is  characterized  by  rich  structure  and 
strong  gradients  in  wind-stress  curl  (Fig.  2(g)  with 
overall  positive  curl  in  the  north-central  and 
eastern  part  of  the  basin  and  negative  curl  along 
the  western  side  of  the  basin,  except  that  the  latter 
is  bisected  by  a  corridor  of  positive  curl  near  the 
East  Korea  Bight. 

Six-hourly  surface  wind-stresses  over  1979-95 
from  a  reanalysis  performed  with  the  National 
Centers  for  Environmental  Prediction  (NCEP) 
atmospheric  model  and  data  assimilation  system 
were  used  to  create  the  final  monthly  climatology 
used  here.  It  is  an  improvement  over  the  NCEP 
operational  nowcast  product  for  reasons  similar  to 
the  case  of  ECMWF,  including  more  atmospheric 
observations  (Kalnay  et  a!.,  1996).  The  analysis 
was  computed  at  spectral  T62  and  subsequently 
transformed  to  a  2.5  x  2.5  grid.  In  the  NCEP  re¬ 
analysis,  the  surface  wind-stress  was  computed 
using  a  drag  coefficient  that  is  a  function  of  the 
bulk  Richardson  number  (Townsend  et  ah,  2000). 
Although  this  data  set  has  relatively  coarse 
resolution,  the  wind-stress  curl  depicts  small-scale 
structure  with  significant  variation  in  magnitude 
(note  from  Fig.  3  that  this  data  set  has  the  second 
highest  basin-averaged  wind  energy).  The  overall 
pattern  is  characterized  by  positive  (negative)  curl 
in  the  eastern  (western)  part  of  the  basin. 

All  of  the  wind-forcing  data  sets  exhibit 
monsoonal  seasonality,  with  weak  winds  from 
the  southeast  in  summer  and  strong  winds  from 
the  northwest  in  winter,  although  significant  local 
variations  in  the  strength  and  pattern  of  wind- 
stress  and  wind-stress  curl  arc  present.  During 
wintertime,  the  mean  synoptic-scale  weather  pat¬ 
terns  are  dominated  by  the  Siberian  high,  which 
produces  north  and  northwest  winds  traveling 
from  the  Asian  continent  over  the  JES,  often  in  the 
form  of  cold-air  outbreaks.  The  seasonality  in  the 
wind-forcing  is  reflected  in  the  seasonal  signal  of 
wind  energy  in  all  of  the  wind-forcing  products 
(Fig.  3).  The  wind  energy  is  computed  as  the 
root  mean  square  of  tv  and  xv.  On  average,  Fig.  3 
shows  that  the  NOGAPS  forcing  has  the  strongest 
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I  ig.  3.  The  iinnual  cycle  of  the  wind  energy  in  N  m~  averaged  per  month  over  the  JES  computational  domain  for  the  seven  different 
wind- forcing  data  sets.  Sec  key  for  explanation. 


basin-wide  wind  energy,  while  Na  has  the  lowest. 
In  general,  this  is  consistent  with  the  magnitude  of 
the  wind-stress  curl  values  (Fig.  2),  although  some 
data  sets  contain  relatively  low  values  of  wind 
energy  (such  as  ECI0M  and  EC  1000),  but  have 
relatively  high  magnitudes  of  curl  and  gradients  of 
curl. 


4.  Linear  response  to  the  wind-forcing  products 

The  linear  solutions  produced  from  a  variety  of 
wind  sets  provide  a  clear  lowest-order  response  to 
the  wind-stress  forcing.  As  such,  these  solutions 
provide  insight  into  the  dynamics  of  observed 
circulation  features  as  well  as  insight  into  non¬ 
linear  model  response  to  the  choice  of  wind¬ 
forcing  product. 


Following  the  approach  of  HH2000,  the  linear 
response  to  each  of  the  atmospheric  forcing  data 
sets  was  calculated  using  a  1.5-layer  reduced 
gravity  model.  These  simulations  represent  the 
lowest-order  dynamics  as  they  have  one  internal 
vertical  mode  and  a  bottom  layer  that  is  infinitely 
deep  and  at  rest.  For  these  simulations,  linearity 
was  achieved  by  reducing  the  amplitude  of  the 
wind  and  straits  forcing  by  a  factor  of  103.  This 
has  the  same  effect  as  excluding  the  nonlinear 
terms  in  the  governing  equations  without  having 
to  change  the  ocean  model  source  code.  HH2000 
examined  the  linear  response  due  to  HR  wind  and 
straits  forcing  both  separately  and  collectively,  but 
here  all  the  simulations  include  the  straits  forcing. 
The  linear  simulations  are  characterized  by  very 
small  local-time  derivatives,  small  nonlinear  ad- 
vective  terms,  and  small-amplitude  variations  in 
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Fig.  4.  Mean  SSH  and  currents  from  the  linear  solutions  forced  by  (A)  throughflovv  only,  throughflow  and  monthly  wind-stress  forcing 
by  (B)  HR,  (C)  COADS,  (D)  NA.  (E)  EC10M.  (F)  ECI000,  (G)  NOGAPS.  (H)  NCEP.  The  solutions  were  linearized  by  reducing  the 
amplitude  of  the  throughflow  and  wind  forcing  by  ICT 


the  interface  depth.  Thus,  the  possibility  of 
diapycnal  mixing  is  excluded. 

Fig.  4  shows  the  annual  mean  SSH  and  currents 
from  eight  1.5-layer,  reduced-gravity  simulations 
forced  with  different  monthly  mean  wmd-forcing 
products,  except  for  Fig.  4a,  which  is  forced  only 
with  flow  through  the  straits  (with  the  amplitude 
reduced  by  103).  Dynamically,  all  of  the  linear 
solutions  in  Fig.  4  are  essentially  the  same  as 
Munk  (1950),  with  a  Sverdrup  (1947)  interior 
calculated  using  realistic  wind  and  throughflow 
forcing  and  coastline  geometry.  However,  in  these 
simulations  horizontal  friction  (A  =  50irr/s)  is 
applied  everywhere  and  the  SSH  field  includes  a 
component  that  balances  the  wind-stress  and 
allows  departures  from  constant  SSH  along  the 
boundary  (even  where  there  is  no  flow  through  the 


boundary),  but  this  component  does  not  contri¬ 
bute  to  the  velocity  field. 

Fig.  4a  depicts  the  linear  response  to  the 
throughflow  forcing  only,  and  as  such  serves  as  a 
baseline  simulation  for  comparison  to  those  that 
include  the  different  wind-forcing  data  sets. 
Although  it  is  uncommon  to  consider  a  model 
result  driven  by  open-boundary  conditions  (only) 
in  the  context  of  Sverdrup  dynamics,  those 
dynamics  are  still  valid,  and  in  this  case  the 
TWC  is  manifested  as  a  Munk  (p1'3)  boundary 
layer  that  hugs  the  western  boundary  until  it  is 
forced  to  separate  from  the  coast  of  Korea  and 
flow  eastward  to  satisfy  the  constraints  of  flow  out 
of  the  Tsugaru  and  Soya  Straits.  That  is  why  these 
zonal  flows  separate  from  the  coast  at  the  latitude 
of  the  straits,  with  two-thirds  of  the  flow  exiting 
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through  Tsugaru  and  the  remaining  one-third 
through  Soya  (a  distribution  that  is  consistent 
with  the  open  boundary  configuration  described  in 
Section  2).  The  grossest  deficiency  in  this  solution 
is  the  northward  flow  along  the  Siberian  coast 
(opposite  to  what  is  observed)  and  the  lack  of  any 
How  to  the  northeast  along  the  coast  of  Honshu. 
As  described  in  HH2000,  the  resulting  SSH 
patterns  from  the  simulations  forced  with  different 
wind  climatologies  are  a  linear  superposition  of 
the  throughflow-driven  solution  (Fig.  4a)  and 
the  solution  produced  by  wind  forcing  only  (see 
Fig.  3  in  HH2000). 

The  mean  SSH  shown  in  Fig.  4b  h  is  the  linear 
solution  for  the  simulations  forced  with  HR, 
GOADS,  Na,  EClOm,  EC1000,  NOGAPS,  and 
NCEP,  respectively,  where  the  abbreviations  are 
defined  in  Section  3.  The  most  prominent  differ¬ 
ence  between  Fig.  4a  and  Fig.  4b-h  is  the  presence 
of  cyclonic  circulation  north  of  the  subpolar  front, 
consistent  with  observations.  This  suggests  that 
the  LCC  is  a  predominantly  wind-driven  current 
that  can  exist  in  the  absence  of  buoyancy  effects. 
All  show  a  Munk  (p1  3)  boundary  layer  that  hugs 
the  western  boundary  east  of  Korea  except  for  the 
case  forced  with  NCEP,  which  produces  the  most 
unrealistic  flow  patterns  of  all  the  wind  sets.  All 
simulations  that  include  wind  forcing  except  the 
NCEP-forccd  simulation  also  show  an  unrealistic 
south  west  ward  return  flow  along  at  least  part  of 
the  Honshu  north  coast. 


5.  Nonlinear  response  to  the  wind-forcing  products 

The  linear  simulations  forced  with  winds  and 
through  How  discussed  in  the  previous  section 
demonstrated  that  many  of  the  large-scale  circula¬ 
tion  features  in  the  JES  could  be  grossly  repro¬ 
duced  with  minimal  dynamical  ingredients  when 
the  model  is  forced  with  winds  and  throughflow. 
In  particular,  all  simulated  the  cyclonic  circulation 
of  the  subpolar  gyre,  but  most  had  unrealistic 
southwestward  flow  along  the  coast  of  Honshu.  In 
the  following  section,  the  model  sensitivity  to  the 
different  wind-forcing  data  sets  is  examined  in 
simulations  with  more  complexity.  The  changes 
include  nonlinearity,  the  addition  of  bottom 


topography,  multiple  internal  modes,  diapycnal 
mixing,  and  ventilation  of  layer  interfaces.  These 
additions  allow  the  possibility  of  barotropic/ 
baroclinic  flow  instabilities,  the  formation  of 
inertial  boundary  layers  and  jets,  and  isopycnal 
outcropping. 

Fig.  5a  shows  the  mean  SSH  for  the  simula¬ 
tion  with  straits  throughflow  forcing  only,  and 
Figs.  5b-h  show  the  mean  SSH  from  the  simula¬ 
tions  forced  with  the  seven  different  monthly 
wind-stress  climatologies  in  addition  to  straits 
throughflow.  The  case  forced  with  straits  through- 
flow  only  (Fig.  5a)  is  able  to  simulate  almost  all  of 
the  major  current  systems  remarkably  well,  the 
exception(s)  being  the  NKCC  and  the  unrealistic 
“overshoot”  of  the  EKWC.  Although  some 
observations  (Shin  et  al.,  2001;  Suh  et  al.,  1999) 
suggest  the  presence  of  a  warm  eddy  in  the  East 
Korea  Bay,  it  is  unclear  if  this  eddy  is  formed 
locally  by  wind  stress  curl,  or  as  a  pinched-off  eddy 
of  the  EKWC,  or  by  trapped  cold  water  sur¬ 
rounded  by  a  migrating  warm  filament  associated 
with  the  warmer  EKWC  water,  or  some  other 
mechanism.  Tn  this  paper,  overshoot  of  the 
EKWC  is  when  the  mean  EKWC  flows  northward 
continuously  along  the  Korean  coast  past  ^40  N. 
It  is  noteworthy  that  the  straits-forced  only 
simulation  reproduces  two  major  circulation  fea¬ 
tures  not  seen  in  the  linear  simulation  (Fig.  4a). 
One  is  the  southward-flowing  LCC  and  the  overall 
cyclonic  flow  of  the  subpolar  gyre  (albeit  weakly) 
in  the  absence  of  any  wind  forcing,  as  these 
features  are  believed  to  be  connected  with 
the  wind-driven  circulation  in  the  JES.  In  this 
case,  they  are  associated  with  a  weak  nonlinear 
recirculation  gyre.  The  second  feature  is  the 
northeastward-flowing  NB  of  the  TWC  along  the 
north  coast  of  Honshu.  In  this  simulation,  the  NB 
is  driven  entirely  by  isopycnal  outcropping  as 
discussed  in  HH2000.  This  is  the  case  for  all  of  the 
nonlinear  simulations  in  Fig.  5,  where  the  corre¬ 
sponding  linear  simulation  shows  no  wind-driven 
contribution  to  the  NB,  i.e.  the  simulations  forced 
by  throughflow  only  and  throughflow  plus  HR, 
COADS^Na,  and  EC  1000. 

The  dynamics  of  the  linear  and  nonlinear  1/8 
simulations  forced  by  throughflow  only  and 
throughflow  plus  HR  are  discussed  extensively  in 
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Fig.  5.  Mean  SSH  (color)  and  enrrenls  (vectors)  from  four-layer  1 '8  nonlinear  simulations  forced  by  (A)  throughflow  only, 
ill  rough  flow  and  monthly  wind-stress  forcing  by  (B)  HR.  (C)  COADS.  (D)  NA.  (H)  liCI0M,(F)  ECI000,  (G)  NOGAPS.  (H)  NCEP. 
Contour  interval  is  1  cm. 


HH2000.  With  these  forcings,  both  the  linear  and 
nonlinear  1/8  simulations  show  overshoot  of  the 
EKWC  separation  from  the  Korean  Coast. 
HH2000  also  discuss  the  1/8  nonlinear  simulation 
forced  by  throughflow  plus  EC10M  which  shows 
realistic  EKWC  separation,  but  not  a  correspond¬ 
ing  linear  simulation.  Thus,  they  did  not  consider 
the  substantial  role  that  the  NB  may  play  in 
EKWC  separation,  a  topic  discussed  shortly. 

Figs.  5b  and  c  show  the  mean  SSH  from  the 
simulations  forced  with  the  HR  and  COADS 
monthly  wind  climatologies,  respectively.  In  these 
simulations  the  mean  SSH  is  quite  similar  to  that 
shown  in  Fig.  5a,  except  that  the  strength  of  the 
subpolar  gyre  is  enhanced  because  of  positive 
wind-stress  curl  in  that  region.  However,  the 


EKWC  still  overshoots  the  observed  separation 
latitude  at  1/8  resolution  because  of  negative  or 
relatively  weak  positive  wind-stress  curl  in  the 
western  subpolar  gyre  region  (Fig.  2b-c).  In 
contrast,  the  mean  SSH  from  the  simulation  with 
atmospheric  forcing  from  the  Na  monthly  wind- 
stress  climatology  is  shown  in  Fig.  5d.  In  this 
simulation,  the  LCC  continues  southward  as  the 
NKCC  until  its  confluence  with  the  EKWC  at 
about  38  N,  consistent  with  the  observed  separa¬ 
tion  latitude  of  the  EKWC.  Also  consistent  with 
the  observations,  part  of  the  EKWC  flow  con¬ 
tinues  eastward  along  the  subpolar  front,  but  some 
recirculates  to  the  south  before  flowing  north¬ 
eastward  along  the  coast  of  Honshu.  Models 
forced  with  this  wind  set  (Seung,  1992;  Seung 
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and  Kim,  1993;  Kim  and  Yoon,  1996,  1999), 
including  those  with  relatively  coarse  horizontal 
grid  resolution,  tend  to  show  the  EKWC  leaving 
the  coast  near  the  observed  separation  latitude 
(although  some  may  include  additional  dynamics 
or  forcing  terms  such  as  buoyancy  forcing).  This  is 
due  to  a  local  maximum  in  the  amplitude  of 
positive  wind-stress  curl  near  East  Korea  Bay  and 
the  large  area  of  positive  wind-stress  curl  in  the 
northern  region  of  the  basin  (Fig.  2d),  even  though 
the  total  wind  energy  is  the  lowest  of  all  the  wind¬ 
forcing  data  sets  examined  (Fig.  3).  In  the  linear 
simulation  with  Na  wind  forcing,  note  that  the 
EKWC  bifurcates  into  a  wind  forced  branch  that 
leaves  the  coast  between  38r  and  39  N  and  a 
through  flow  forced  branch  that  overshoots  and 
separates  near  41  N  (Fig.  4d).  In  the  nonlinear 
solution,  part  of  the  flow  required  to  feed  the 
outflow  straits  is  diverted  into  the  NB,  making 
it  easier  for  the  wind-stress  curl  to  dominate 
in  determining  the  EKWC  separation  latitude 
(Fig.  5d).  A  similar  situation  exists  for  the 
ECIOM  linear  and  nonlinear  simulations. 

Figs.  5c  and  f  show  the  mean  SSH  from  the 
nonlinear  simulations  forced  with  the  monthly 
climatologies  formed  from  the  ECMWF  re-analy¬ 
sis  (10  m  and  1000  mb,  respectively).  Recall  that 
the  ECIOM  winds  have  1.125°  resolution  and 
EC  1000  have  2.5  resolution.  Here  we  assume  the 
primary  differences  in  the  structure,  magnitude 
and  gradients  of  the  wind-stress  curl  in  the  EC  1000 
data  set  arc  due  to  truncation  of  the  spectral 
grid  to  T47  (and  subsequent  interpolation  to  a 
2.5  x  2.5  Mercator  grid).  A  notable  feature  in  the 
ECIOM  data  set  is  the  zonal  corridor  of  positive 
wind-stress  curl  extending  eastward  from  East 
Korea  Bay  at  39  N  (which  is  also  seen  in 
the  NOGAPS  and  Na  wind-stress  curl  patterns 
(  Figs.  2d  e.g)).  As  a  result,  the  mean  separation 
latitude  of  the  EKWC  is  farther  to  the  south  in  the 
simulation  forced  with  ECIOM  (Fig.  5e)  than  the 
N-a- forced  simulation. 

Another  significant  feature  of  the  ECIOM  mean 
wind-stress  curl  is  the  strong  region  of  negative 
wind -stress  curl  in  the  central  JES  centered  near 
40  N,  134  E  Hanked  by  regions  of  strong  positive 
wind-stress  curl  with  a  band  of  strong  negative 
wind-stress  curl  along  the  eastern  boundary  in  the 


corresponding  latitude  range.  This  results  in  a 
pattern  of  large  meanders  across  the  basin  between 
38°N  and  40°N  in  the  linear  simulation  (Fig.  4e). 
These  are  not  seen  in  any  of  the  other  simulations, 
including  EC  1000.  A  similar  pattern  of  meanders 
is  seen  in  the  nonlinear  simulation,  but  with  the 
northward  meander  in  the  central  basin  shifted 
eastward  2  so  that  it  lies  over  the  Yamato 
Rise.  This  eastward  shift  is  due  to  the  mean 
abyssal  flow  along  the  geostrophic  contours  of  the 
topography.  The  role  of  abyssal  currents  in 
steering  upper  currents  is  discussed  in  Section  6. 
Another  significant  difference  between  the  linear 
and  nonlinear  simulations  with  ECIOM  wind 
forcing  is  seen  in  the  pathway  of  the  subpolar 
front,  a  dual  pathway  in  the  linear  simulation  (the 
meandering  pathway  and  a  more  zonal  pathway 
along  41  N)  and  only  the  meandering  pathway  in 
the  nonlinear  simulation.  This  is  a  consequence  of 
the  dynamics  used  in  explaining  EKWC  separa¬ 
tion  from  the  coast  in  the  nonlinear  Na  simulation 
(Fig.  5d),  where  much  of  the  flow  connecting  the 
inflow  and  outflow  ports  is  shifted  to  the  NB. 

In  marked  contrast  to  ECIOM,  the  EC  1000 
wind-stress  curl  is  very  smooth  and  dominated  by 
positive  wind-stress  curl  in  the  northeast  half  of 
the  basin  and  negative  in  the  southwestern  half. 
Gone  is  the  region  of  strong  positive  wind-stress 
curl  east  of  East  Korean  Bay  that  played  a  key  role 
in  EKWC  separation,  but  a  remnant  of  positive 
wind-stress  curl  remains  in  East  Korea  Bay.  Also 
gone  is  most  of  the  meandering  in  the  subpolar 
frontal  region  of  the  linear  and  nonlinear  simula¬ 
tions  with  EC1000  wind  forcing,  including  the 
northward  meander  and  gyre  above  the  Yamato 
Rise.  However,  the  anticyclonic  gyre  east  of  Korea 
is  the  strongest  of  all  of  the  simulations.  The  linear 
simulation  shows  a  single  EKWC  pathway  that 
leaves  the  coast  between  40  N  and  41  N.  There  is 
no  bifurcation  into  the  wind-driven  and  through- 
flow-driven  pathways  as  seen  with  Na  and  ECIOM 
wind  forcing.  Yet,  the  nonlinear  simulation  shows 
partial  separation  of  the  EKWC  starting  at  38°N 
and  complete  separation  by  38. 8:  N.  Here  we 
investigate  the  hypothesis  that  this  separation 
is  inertial  and  can  be  described  by  a  constant 
absolute  vorticity  (CAV)  trajectory  at  least  to 
the  northernmost  latitude  of  the  mean  EKWC 
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trajectory  after  separation.  In  addition,  the  posi¬ 
tive  wind-stress  curl  in  East  Korea  Bay  and 
anticyclonic  abyssal  flow  around  the  rise  near  the 
separation  point  may  have  some  influence  on  the 
pathway  after  separation.  HH2000  found  that 
flow  around  this  rise  could  play  a  key  role  in 
EKWC  separation  from  the  coast.  Of  the  simula¬ 
tions  in  Fig.  5,  the  simulation  using  EC1000  wind 
forcing  has  the  strongest  flow  around  this  rise  with 
a  maximum  greater  than  I  cm/s  near  129.5  E, 
38.2  N,  but  still  much  weaker  than  in  simulations 
where  this  had  an  effect  in  HH2000. 

A  CAV  trajectory  occurs  in  a  special  case  where 
conservation  of  potential  vorticity  reduces  to 
CAV.  It  occurs  at  the  core  of  a  steady,  geos- 
trophic,  inviscid  current  (no  vorticity  due  to 
shear).  In  a  1.5-layer  nonlinear  reduced-gravity 
model  on  a  fi-plane,  it  forms  a  streamline  with 
constant  speed  along  a  constant  layer  thickness 
contour.  Starting  at  an  inflection  point  in  the 
trajectory,  the  maximum  northward  penetration  of 
a  CAV  trajectory  is  given  by 

b  =  cl  sin  !©(),  ( 1 ) 

where 

d  =  2(Vc/fi)''2  sin  l/20o  (2) 

and  Vc  is  the  speed  at  the  core  of  the  current, 
P  =  df  /d  v\  and  0O  is  the  angle  of  the  flow  at  the 
inflection  point  measured  counterclockwise  from 
due  east  (Hurlburt  and  Thompson,  1982).  Of  the 
eight  JES  simulations  shown  in  Fig.  5,  the 
simulations  with  EC10M  and  EC  1 000  wind  for¬ 
cing  have  the  strongest  EKWC  at  the  separation 
point  with  Fc  =  36cm/s  in  the  top  layer  and  are 
the  two  where  CAV  trajectory  theory  may  be  most 
useful  and  applicable.  However,  at  the  northern¬ 
most  point  of  the  trajectory,  Vc  drops  to  23  cm/s  in 
the  simulation  using  ECIOM  and  25 cm/s  in  the 
simulation  using  EC  1000.  Hence,  we  choose  to 
specify  the  other  parameters  and  solve  for 

[/c  =  /?/72/4^sin  ie0)  .  (3) 

It  is  also  difficult  to  precisely  determine  the 
latitude  of  the  inflection  point,  a  potential  source 
of  significant  error,  since  it  is  used  in  calculating 


ft2.  We  used  the  separation  latitude  (37.25  N), 
P=  1.77  x  10“"  m-'s-'  at  38'N,  and  0O  = 
104.6  (106.5  )  to  obtain  ft  =  2.33  km  (272  km) 

from  the  simulations  and  Vc  =  38cm/s  (51  cm/s) 
from  Eq.  (3)  for  the  simulations  using  ECIOM 
(EC1000).  The  values  of  Vc  obtained  from  CAV 
trajectory  theory  are  obviously  higher  than  in  the 
numerical  simulations,  particularly  for  the  simula¬ 
tion  using  EC  1 000.  Continuation  of  some  fric¬ 
tional  contact  of  the  current  with  the  boundary 
layer  after  separation,  more  strongly  in  the 
simulation  using  EC1000,  is  an  explanation  that 
is  consistent  with  these  results.  Still,  the  results 
indicate  that  the  CAV  dynamics  are  a  useful 
benchmark  in  assessing  the  deviation  from  their 
prediction,  as  will  be  seen  in  Section  6. 

NOGAPS  and  NCEP  are  the  strongest  wind- 
forcing  products  and  tend  to  dominate  over 
the  throughflow  forcing  in  the  model  simulations 
(Fig.  2).  The  NCEP  pattern  of  mean  wind-stress 
curl  is  largely  the  opposite  of  Na  (and  much 
stronger)  with  negative  wind-stress  curl  over  the 
northwestern  half  of  the  basin  and  positive  in 
the  southeastern  half.  The  simulations  forced  with 
the  NOGAPS  wind  stresses  (Fig.  5g)  and  the 
NCEP  reanalysis  (Fig.  5h)  are  characterized  by 
strong  cyclonic  flow  over  a  large  part  of  the  JES 
basin.  The  pattern  and  magnitude  of  the  wind- 
stress  curl  in  these  two  data  sets  result  in  strong 
Ekman  suction  in  the  subpolar  gyre,  enhancing  its 
strength.  In  fact,  in  both  of  these  simulations  the 
subpolar  front  exists  only  to  the  extent  that  it 
forms  the  southern  boundary  of  the  subpolar  gyre. 
Accordingly,  the  NB  is  stronger  in  these  two 
solutions  than  in  any  of  the  others.  Indeed,  one 
could  argue  that  these  simulations  reproduce  the 
separation  of  the  EKWC  from  the  Korean  coast 
and  northward  flow  along  the  coast  of  Honshu 
reasonably  well.  A  more  detailed  discussion  of  the 
EKWC  and  the  current  pathway  after  separation 
from  the  coast,  in  relation  to  the  1/32'  simulations 
using  these  two  wind-forcing  products,  is  con¬ 
tained  in  Section  6.  However,  other  aspects  of 
these  solutions,  particularly  the  double  gyre 
cyclonic  circulation  feature  near  133  E,  38  N  and 
the  strong  anticyclonic  eddy  near  the  East  Korea 
Bay  in  Fig.  5h  are  less  realistic.  These  two  features 
are  due  to  the  region  of  strong  positive  wind-stress 
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curl  in  the  southernmost  JES  and  strong  negative 
wind-stress  curl  cast  of  East  Korea  Bay,  respec¬ 
tively.  in  the  NCEP  forcing  (Fig.  2h). 

Although  the  wind  stress  used  by  NLOM  acts 
only  on  the  layer  in  direct  contact  with  the 
atmosphere,  the  mean  abyssal  circulation  is  also 
affected  by  the  surface  wind  forcing,  in  that  wind 
forcing  data  sets  with  stronger  wind-stress  and 
positive  wind-stress  curl  (e.g.  Figs.  2  and  3)  tend  to 
promote  more  isopycnal  outcropping  and  thus 
more  vigorous  diapycnal  mixing  and  exchange 
between  the  surface  and  subsurface  layers,  includ¬ 
ing  the  abyssal  layer.  This  is  reflected  in  maps  of 
abyssal -layer  kinetic  energy  of  the  mean  flow 
(KEM)  as  shown  in  Fig.  6.  As  such,  the  simulation 
forced  with  NOGAPS,  which  has  the  strongest 


positive  wind-stress  curl  with  a  relatively  large 
zonal  extent,  has  a  relatively  strong  cyclonic  mean 
abyssal  circulation  (Fig.  6g).  This  is  consistent 
with  limited  but  recent  observations,  which  suggest 
deep  cyclonic  flow  in  the  JES,  but  with  more 
episodic  events  (rather  than  mean)  of  the  magni¬ 
tude  seen  in  Fig.  6g  (Takematsu  et  al.,  1999; 
Senjyu  et  ah,  2002).  Conversely,  the  cases  forced 
with  EC1000,  EClOm,  and  NCEP  show  weaker 
cyclonic  abyssal  circulation,  and  the  cases  forced 
with  HR,  COADS,  and  inflow/outflow  only  have 
weakly  anticyclonic  abyssal  circulation.  The  rela¬ 
tively  strong  mean  abyssal  circulation  in  the  1/8  - 
NOGAPS-forced  simulation  is  driven  by  relatively 
strong  vertical  mixing  in  the  subpolar  gyre,  where 
there  is  isopycnal  outcropping.  This  is  different 


Fig.  6.  Kinetic  energy  per  unit  mass  of  the  mean  abyssal  flow  (color)  and  mean  abyssal  currents  (vectors)  for  the  1/8  simulations 
shown  in  Fig.  5:  (A)  ihroughflow  only,  ihroughflow  and  monthly  wind-stress  forcing  by  (B)  HR.  (C)  COADS.  (D)  NA,  (E)  ECI0M, 
(I  )  LCIOOO,  (G)  NOGAPS,  (H)  NCEP 
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from  the  mixed  barotropic/baroclinic  instability 
mechanism  described  in  HH2000  that  energizes  the 
abyssal  layer,  but  the  result  is  basically  the  same. 
Both  depict  robust  abyssal-layer  currents  that 
are  topographically  constrained  and  influence 
the  pathways  of  surface-layer  currents  via  upper 
ocean-topographical  coupling,  including  an  influ¬ 
ence  on  the  separation  latitude  of  the  EKWC  in 
some  cases. 

The  location  of  diapycnal  mixing  between  the 
surface  and  abyssal  layers  helps  to  elucidate  the 
location  and  process  of  deep-water  formation  and 
abyssal-layer  current  pathways.  As  described  in 
HH2000,  in  NLOM  outcropped  layers  do  not 
attain  zero-layer  thickness,  but  rather  outcrop  to  a 
specified  layer  thickness  that  controls  activation  of 
vertical  mixing.  The  location  where  the  laycr(s) 
outcrop  is  co-located  with  the  location  of  the 


interfaces  ventilating.  HH2000  describe  the  venti¬ 
lation  patterns  for  the  1/8C  case  forced  with  the 
HR  monthly  wind  climatology,  inflow/outflow 
through  the  straits  (no  wind  forcing),  and  both 
as  a  function  of  the  fraction  of  time  that  a  given 
interface  ventilates.  Those  calculations  have  been 
extended  to  the  simulations  discussed  here  for  the 
deepest  (third)  interface  (Fig.  7).  In  all  simulations 
the  deepest  interface  ventilates  north  of  the 
subpolar  front  but  the  areal  extent  varies  signifi¬ 
cantly,  and  in  all  cases  the  location  where  the 
bottom  interface  is  ventilating  is  largely  co-located 
with  the  pattern  of  positive  wind-stress  curl.  The 
positive  wind-stress  curl  invokes  vigorous  Ekman 
suction,  causing  the  isopycnals  to  dome  upward 
and  outcrop  once  the  specified  minimum  layer 
thickness  is  reached.  Likely  regions  of  intermediate 
and  deep-water  formation  are  located  where  the 


Fig.  7  Maps  of  ihe  percent  of  lime  that  the  third  (bottom)  interface  ventilates  for  the  l  8  throughflow  only  and  seven  different  wind 
climatology  (including  throughflow)  cases.  Units  are  percent  of  time:  contour  interval  is  10%.  Note  that  0%  is  a  separate  color 
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bottom  interface  is  ventilating  most  of  the  time,  In 
all  eases  this  is  north  of  the  subpolar  front  in  the 
subpolar  gyre.  This  is  reasonable,  and  consistent 
with  limited  observations  (e.g.,  Senjyu  and  Sudo, 
1994;  Scung  and  Yoon,  1995),  but  there  are 
substantial  differences  in  the  outcropping  patterns 
depicted  in  Fig.  7.  Hence,  when  using  a  numerical 
model  to  investigate  water-mass  formation  and 
movement,  particular  care  should  be  given  to  the 
choice  of  atmospheric  forcing. 

6.  Comparison  of  nonlinear  model  responses  to  the 
combined  effects  of  wind-forcing  choice  and  upper 
ocean-topographical  coupling 

As  previously  mentioned,  the  simulations  de¬ 
scribed  in  Section  5  were  run  with  1/8'  horizontal 
grid  resolution  to  elucidate  more  clearly  the  impact 
of  wind  forcing  rather  than  flow  instabilities.  This 
is  not  to  say  those  simulations  arc  devoid  of  flow 
instabilities,  as  they  arc  eddy-resolving.  However, 
HH2000  demonstrated  that  the  flow  instabilities 
arc  much  more  pervasive  at  1/32°  resolution,  and 
that  flow  instabilities  at  that  resolution  have  a 
profound  impact  on  the  surface  and  abyssal 
circulation  when  the  model  is  forced  with  HR 
climatological  winds  and  throughflow  forcing.  In 
particular,  the  stronger  baroclinic  instability 
transfers  much  greater  eddy  kinetic  energy  to  the 
abyssal  layer.  This  results  in  much  stronger  abyssal 
flow  and  much  stronger  upper  ocean-  topogra¬ 
phical  coupling  via  a  mixed  baroclinic/barotropic 
instability,  where  baroclinic  instability  is  key  to  the 
coupling  because  of  the  strong  downward  energy 
transfer. 

In  this  section  we  compare  the  mean  surface 
(Fig.  8)  and  mean  abyssal  (Fig.  9)  results  from 
1/32  J  ES  simulations  forced  by  six  of  the  wind 
products.  Simulations  forced  by  COADS  and 
throughflow  only  are  omitted  due  to  similarity  of 
the  1/8  nonlinear  simulations  with  the  1/8  HR- 
forced  simulation.  The  six  simulations  are  identi¬ 
cally  configured  except  for  the  wind  forcing, 
and  all  include  the  throughflow.  They  are  also 
configured  identically  to  the  corresponding  l/8~ 
simulations,  except  for  (1)  the  grid  resolution 
(and  corresponding  improvements  in  resolving  the 


model  boundary  and  bottom  topography)  and  (2) 
the  eddy  viscosity,  which  was  reduced  from  A  = 
50m2/s  at  1/8°  to  5m2/s  at  1/32  . 

Substantial  impacts  are  seen  in  both  the  surface 
layer  (Fig.  5  vs.  Fig.  8)  and  the  abyssal  layer  (Fig.  6 
vs.  Fig  9)  due  to  the  resolution  increase  (and 
associated  eddy  viscosity  decrease)  and  the  strong 
upper  ocean-  topographical  coupling  via  baroclinic 
instability.  Although  the  mean  SSH  of  the  simula¬ 
tions  with  the  different  wind  forcings  generally 
demonstrate  increased  similarity  at  l/32c  resolu¬ 
tion,  the  choice  of  wind-forcing  product  retains  a 
large  impact.  The  mean  surface-layer  circulation 
from  the  l/32r  NOGAPS-forced  simulation  shows 
the  least  change  from  its  1/8  counterpart  because 
the  1/8°  simulation  already  had  relatively  strong 
abyssal  flow  that  steered  upper-ocean  currents, 
including  the  pathway  of  the  subpolar  front  in  the 
vicinity  of  the  Yamato  Rise.  That  abyssal  flow  was 
due  to  greater  isopycnal  outcropping  than  in  the 
other  1/8°  simulations  (as  discussed  in  Section  5). 
However,  the  changes  that  did  occur  in  the  1/32' 
NOGAPS-forced  simulation  improved  it  to  the 
point  that,  unlike  at  1/8  resolution,  it  is  no  longer 
an  outlier  in  relation  to  the  other  1/32'  simulations 
(which  have  also  converged  toward  the  1/32' 
NOGAPS-forced  simulation).  Only  the  NCEP- 
forced  simulation  remains  a  gross  outlier  at  1/32 
resolution. 

The  1/8°  and  1/32'  EC10M  simulations  also 
have  relatively  similar  mean  surface-layer  circula¬ 
tion  because  the  flow  around  the  Yamato  Rise  is 
supported  by  a  feature  in  the  wind-stress  curl  in 
addition  to  influence  by  the  Yamato  Rise  (as 
discussed  in  Section  5).  All  the  1/32  simulations, 
except  for  the  NCEP-forced,  show  strong  impact 
of  the  Yamato  Rise  on  the  pathway  of  the 
subpolar  front,  but,  except  for  the  EClOm 
simulation,  this  is  entirely  due  to  abyssal  current 
steering  of  the  surface-layer  mean  current.  At  1/8 
these  abyssal  currents  were  driven  by  vertical 
mixing  due  to  isopycnal  outcropping,  but  at  1/32 
they  are  driven  almost  entirely  by  downward 
transfer  of  energy  due  to  baroclinic  instability, 
except  for  a  significant  contribution  from  vertical 
mixing  in  the  NOGAPS-forced  simulation  (Fig.  6 
vs.  Fig.  9).  The  l/32;  NCEP-forced  simulation 
does  not  show  a  Yamato  Rise  impact  on  the 
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Pig.  8.  Mean  SSH  (color)  and  currents  (vectors)  from  four-layer  132  nonlinear  simulations  forced  by  (A)  throughflow  only, 
through  flow  and  monthly  wind-stress  forcing  by  <B)  HR.  (C)  COADS.  (D)  NA,  (b)  ECIOM.  (F)  EC1000.  (G)  NOGAPS.  (H)  NCEP 
Contour  interval  is  1  cm. 


subpolar  frontal  pathway  because  this  front 
approaches  the  rise  from  the  northwest  and  lacks 
a  frontal  segment  between  the  EKWC  and  the  west 
side  of  the  Yamato  Rise.  The  reasons  for  this  are 
discussed  shortly. 

HH2000  showed  that  upper  ocean -topographi¬ 
cal  coupling  was  essential  for  realistic  simulation 
of  the  EKWC  separation  from  the  coast  (when 


forced  with  HR  winds)  and  to  prevent  the  over¬ 
shoot  of  the  observed  separation  latitude  seen  in 
the  corresponding  1/8  simulation.  They  also 
demonstrated  that  anticyclonic  abyssal  flow 
around  a  rise  centered  at  38.3  N,  129.8  E  was 
crucial  to  the  separation.  In  addition,  there  is 
strong  abyssal  flow  cyclonically  around  the  rim  of 
the  abyssal  plain  at  ^2500-3000  m  depth  (seen  in 
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I  :ig.  9.  Kinetic  energy  per  unit  mass  of  the  mean  abyssal  flow  (color)  and  mean  abyssal  currents  (vectors)  for  the  1/32  simulations 
shown  in  lig.  8:  (A)  throughfiow  only,  ihroughflow  and  monthly  wind-stress  forcing  by  (B)  HR,  (C)  COAOS.  (D)  NA.  (E)  EC10M, 
(!•')  EC 1 000,  (G)  NOGAPS,  (H)  NCEP. 


all  the  1/32  simulations).  This  flow  proceeds 
southward  to  a  point  immediately  east  of  the 
rise  that  is  crucial  to  EKWC  separation  and 
then  abruptly  changes  to  eastward  flow  along 
38.5  38.8  N  toward  the  Yamato  Rise.  After 
EKWC  separation  from  the  coast,  the  surface 
current  proceeds  eastward  as  the  subpolar  front 
and  parallels  this  eastward  abyssal  current  to  the 
Yamato  Rise.  The  same  is  seen  in  the  l/32c 
EClOOO-forced  simulation  which  has  a  surface 
flow  pattern  that  is  very  similar  to  the  l/32  J  HR- 
forecd  simulation,  but  with  generally  larger 
amplitude. 


The  EC10M-,  Na-,  and  NOGAPS-forccd  simu¬ 
lations  all  had  realistic  EKWC  separation  latitudes 
at  l/8f  resolution  due  to  the  wind  forcing,  and 
that  remains  the  case  at  1/32  resolution.  All 
three  have  post-separation  pathways  that  pass 
southeast  of  the  rise  at  38.3' N,  129.8  E,  and  flow 
around  this  rise  is  less  energized  in  these  simula¬ 
tions.  However,  the  northward  penetration  of 
all  three  pathways  is  limited  to  a  lower  latitude 
than  in  the  HR-  and  EClOOO-forced  simulations 
because  they  approach  the  eastward  abyssal 
current  from  the  south  and  because  of  anticyclonic 
abyssal  flow  around  a  rise  centered  at  38.0  N, 
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1 33.3  JB.  In  the  Na-forced  simulation,  the  eastward 
portion  of  the  surface-current  pathway  lies  along 
the  southern  side  of  the  eastward  abyssal  current 
and  north  of  the  rise. 

In  the  1/8C  and  1/32  NOGAPS-  and  1/32 
EClOM-forced  simulations,  the  post-separation 
pathway  forms  a  meander  where  the  pathway 
returns  toward  the  south  after  reaching  maximum 
northward  penetration.  That  pathway  is  nearly  the 
same  for  the  1/8  and  1/32'  NOGAPS-forced 
simulations.  The  EClOM-forced  simulation,  which 
equaled  the  EClOOO-forced  simulation  in  having 
the  strongest  EKWC  at  1/8  resolution,  has 
formed  a  circular  nonlinear  recirculation  gyre 
at  l/32r  resolution  with  a  diameter  of  200  km 
between  speed  maxima  and  with  maximum  mean 
speed  in  the  EKWC  increasing  from  36  to 
>50cm/s  (the  strongest  of  the  1/32  simulations). 
Other  possible  influences  for  change  in  the  EKWC 
separation  at  1/32  resolution  are  (1)  a  southward 
abyssal  current  adjacent  to  the  western  boundary, 
(2)  a  change  in  the  surface-layer  outcropp¬ 
ing  causing  a  change  in  the  strength  of  the 
NB  (discussed  in  Section  5),  or  a  change  in  the 
subpolar  gyre  in  the  vicinity  of  separation,  and  (3) 
the  strength  of  the  EKWC.  The  maximum  mean 
surface-layer  speed  of  the  EKWC  in  cm/s 
are  for  HR  =  40,  EC10M  =  52,  EC1000  =  43, 
Na  -  24,  NOGAPS  =  29  and  NCEP  =  28.  The 
results  indicate  that  any  influence  of  (1)  or  (2)  is 
not  large.  The  size/strength  of  the  subpolar  gyre 
does  not  show  substantial  changes  at  1/32  ,  but 
for  some  simulations  the  southward  NKCC 
does  show  increases  in  maximum  mean  surface- 
layer  speed,  ~25%  with  EC10M.  0%  with 
EC  1000,  70%  with  Na  and  20%  with  NOGAPS 
forcing.  At  1/32'  the  corresponding  speeds  of 
the  NKCC  (cm/s)  arc  3  with  HR,  8  with 
EC10M,  6  with  EC1000,  26  with  Na,  and  18  with 
NOGAPS  forcing.  There  is  no  NKCC  with  NCEP 
forcing,  just  a  small  cyclonic  gyre  spanning 
38-39  N.  For  all  six  simulations,  the  southward 
deep  western  boundary  current  shows  obvious 
increases  at  1/32  (Fig.  6  vs.  Fig.  9).  Despite  the 
preceding  changes,  the  EKWC  separation  docs 
not  change  appreciably  for  the  EC10M-,  Na-, 
and  NOGAPS-forced  simulations.  In  addition, 
the  NCEP-forced  simulation  still  overshoots  at 


1/32  and  the  strength  of  the  southward  western 
boundary  currents  are  not  significant  to  the 
explanation  for  EKWC  separation  in  the  HR- 
and  EClOOO-forced  simulations,  as  clearly  demon¬ 
strated  by  HH2000  for  the  HR-forced  simula¬ 
tion.  Also,  there  is  no  relation  between  the 
strength  of  the  EKWC  and  its  latitude  of  separa¬ 
tion  or  its  post-separation  pathway,  except  for 
some  speed-related  pathway  modification  in  the 
EC10M  simulation  (which  had  the  strongest 
EKWC). 

For  the  1/32  ECI0M  and  the  1/8  and  1/32' 
NOGAPS-forced  simulations,  the  pathway  alter 
separation  was  tested  for  agreement  with  a  CAV 
trajectory  using  the  methodology  introduced  in 
Section  5.  For  the  1/32  EC10M  (l/32r,  l/8r 
NOGAPS)-forced  simulations  we  used  a  separa¬ 
tion  latitude  of  36.7  N,  ft  =  1.79  x  lO^'m-'s'1 
at  37  N  and  0o  =  92.  E  (93.1 90  )  to  obtain  b  = 
155  (155,  1 22)  km  from  the  simulations  and 
calculated  Vc  =  21  (20,  13) cm/s  from  Eq.  (3). 
Since  the  simulated  value  of  Vc  is  a  maximum  of  52 
(29,  26) cm/s  and  the  minimum  for  Vc  anywhere  in 
layers  1  or  2  between  EKWC  separation  and  its 
northernmost  latitude  is  33  (21,  1 5)  cm/s,  the 
pathway  predicted  by  CAV  trajectory  dynamics 
would  extend  farther  north  (larger  value  ofb)  than 
simulated  in  these  three  cases,  the  opposite  of  the 
result  obtained  at  1/8  using  EC10M  forcing 
(Section  5). 

The  agreement  with  CAV  trajectory  theory 
for  the  1/32  NOGAPS-forced  simulation  is 
relatively  good,  but  that  may  not  be  due  to 
CAV  trajectory  dynamics.  Baroclinic  instability 
is  strong  in  the  1/32  simulations  and  can 
shorten  the  meander  wavelength  compared  to  a 
CAV  trajectory.  However,  baroclinic  instability  is 
weak  at  1/8  resolution,  making  this  explanation 
inconsistent  with  the  similarity  between  the  1/8' 
and  1/32  pathways  in  the  NOGAPS-forced 
simulations.  Further  examination  of  the  l/32c 
NOGAPS-forced  simulation  (and  the  other  two 
simulations)  reveals  that,  after  EKWC  separation, 
the  current  is  augmented  by  flow  joining  the 
eastward  jet  from  the  north  via  the  NKCC 
and  a  second  current  (with  a  mean  speed  of 
~1 1  cm/s  in  the  1/32  NOGAPS-forced  simula¬ 
tion)  merging  near  the  northernmost  point  of 
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the  mean  meander  after  separation.  The  latter  is 
strongly  barotropie  (~6cm/s  in  the  abyssal  layer). 
In  the  abyssal  layer,  this  is  the  same  current  that 
turns  eastward  and  constrains  the  mean  pathway 
in  the  1/32  HR-,  EC1000-,  and  Na-forced 
simulations.  Here  it  would  limit  the  northward 
meander  penetration  of  the  separated  jet  in  the 
1/32  NOGAPS-  and  EClOM-forced  simulations, 
which  have  the  same  northward  penetration.  In 
the  1/8'  NOGAPS-forccd  simulation,  the  north¬ 
ward  penetration  is  insufficient  to  be  limited  by 
this  abyssal  current,  which  is  present  at  1-  2cm/$  in 
that  simulation. 

Additional  influence  on  the  pathways  of  the 
1/32  simulations  comes  from  stronger  anticyclo- 
nic  abyssal  flow  around  the  rise  centered  near 
38.0  N,  133.3  E  ( <  I  cm/s  in  the  1/8  NOGAPS- 
forccd  simulation)  and,  in  all  three,  southwest- 
ward  abyssal  flow  entering  the  Ulleung  basin 
via  a  deep  channel  on  the  east  side  of  the  rise. 
About  one-half  degree  to  the  southeast  of  this 
deep  channel  compensating  northeastward  flow 
exits  along  the  eastern  side  of  the  channel.  This 
combination  of  abyssal  flow  steers  an  indentation 
in  the  southeastward  surface  jet  in  all  three 
simulations.  It  is  best  seen  in  the  l/32r  NO- 
G  A  PS-forced  simulation  where  it  is  centered 
near  37.6  N.  131.9  E  (Fig.  9e)  (a  much  larger 
indentation  will  be  shown  in  Section  7).  This 
abyssal  steering  also  can  limit  the  northward 
penetration  of  the  surface  current  meander.  Unlike 
the  simulations  forced  by  the  other  wind  sets 
(except  NCEP),  the  NOGAPS  wind-stress  curl 
strongly  forces  the  current  southward  after  EKWC 
separation  in  the  linear  simulation  (Fig.  4g)  and 
therefore  can  also  play  a  major  role  in  limiting 
northward  penetration  of  the  post-separation 
meander  in  the  nonlinear  NOGAPS-forced  simu¬ 
lations. 

Only  the  NCEP- forced  simulation  continues  to 
show  EKWC  overshoot  at  1/32  resolution.  It  does 
form  an  eddy  (in  the  mean)  over  the  seamount  at 
38.3  N,  129.8  E,  but  there  is  sufficient  negative 
wind-stress  curl  support  that  part  of  the  EKWC 
continues  northward  all  the  way  to  Vladivostok  as 
in  the  1/8  linear  and  nonlinear  simulations  (Figs. 
4h  and  5h),  including  an  offshore  northward 
segment  around  38  N. 


7.  Deterministic  vs.  nondeterministic  response  to 
wind  forcing 

A  series  of  six  interannual  simulations  are  used 
to  quantify  the  relative  influence  of  the  determi¬ 
nistic  response  to  the  wind  and  throughflow 
forcing  vs.  the  nondeterministic  responses  due  to 
flow  instabilities  on  the  JES  circulation.  For  this 
purpose,  the  six  simulations  used  l/32f  resolution 
and  were  forced  with  the  six-hourly  ECMW'F  10  m 
reanalysis  winds  over  the  time  frame  1979-93.  The 
1979-93  ECMWF  reanalysis  also  was  used  to 
create  the  climatological  ECMWF  10m  and 
ECMWF  1000 mb  reanalysis  data  sets  described 
in  Section  3.  The  initial  state  for  the  otherwise 
identical  simulations  was  varied  by  starting  each 
realization  1  month  later  from  the  climatologically 
forced  1/32'  EC10M  simulation  after  it  reached 
statistical  equilibrium.  The  atmospheric  forcing 
also  progressively  started  1  month  later  in  1979. 

Fig.  10  shows  the  annual  mean  SSH  formed 
over  1991  from  each  of  the  six  realizations. 
Because  the  six  simulations  were  identical  except 
for  the  initial  state,  the  differences  between  them 
are  due  to  mesoscalc  flow  instabilities,  and,  as 
Fig.  10  demonstrates,  the  flow  instabilities  have  a 
dramatic  impact  on  the  1-year  mean  SSH  and 
current  pathways  and  their  dynamics.  We  will 
briefly  consider  the  post-EKWC  separation  path¬ 
way  as  an  example.  All  realizations  depict  the 
EKWC  separating  near  37  N,  although  there  is 
substantial  seasonal  variation  (seasonal  means 
not  shown).  Case  1  has  a  post-EKWC  separation 
pathway  similar  to  the  mean  SSH  from  the  1/32' 
HR  and  EC1000  simulations  (Fig.  8),  while  the 
others  have  post-separation  current  pathways 
more  like  the  climatologically  forced  1/327 
EC10M  simulation,  but  with  varying  indentation 
amplitudes  on  the  east  side  of  the  meander  and 
varying  strengths  of  the  nonlinear  recirculation 
gyre.  The  significance  and  dynamics  of  the 
indentations  are  discussed  in  Section  6.  There  is 
no  clear  relation  between  the  indentation  ampli¬ 
tude  and  the  northward  penetration  of  the  mean 
meander,  although  case  3  (which  has  the  largest 
indentation)  has  the  second  lowest  penetration 
latitude.  Case  4  has  the  least  northward  penetra¬ 
tion,  the  strongest  nonlinear  recirculation  gyre. 
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Fig.  10.  Mean  SSH  for  1991  from  six  realizations  of  four-layer  1,32  NLOM  forced  with  synoptic  HCMWF  10m  reanalysis  wind- 
stress.  Fach  realization  differs  only  by  the  initial  slate.  Contour  interval  is  1  cm. 


and  the  strongest  EKWC  and  post-separation 
current.  This  could  be  due  to  abyssal  current 
steering  that  produced  the  indentation  in  case  3 
acting  on  a  more  inertial  surface  current  or  due  to 
wavelength  shortening  (compared  to  a  CAV 
trajectory)  from  baroclinic  instability.  All  but 
cases  3  and  4  have  a  northward  penetration 
latitude  that  would  be  limited  by  the  eastward 
abyssal  current  along  38.5-38.8  N,  a  topic  dis¬ 
cussed  in  Section  6. 

There  are  also  significant  variations  in  the  flow 
patterns  in  the  interior  of  the  basin.  Some  cases  (1 
and  4;  Figs.  10a  and  lOd)  depict  a  relatively  strong 


anticyclonic  eddy  north  of  the  subpolar  front,  a 
feature  that  to  date  is  not  supported  by  observa¬ 
tions  (to  the  authors'  best  knowledge).  There  is 
substantial  variability  in  the  position  and  size  of 
the  persistent  (I -year  mean)  eddies,  particularly  in 
the  southern  part  of  the  basin.  The  differences  in 
the  solutions  shown  in  Fig.  10  are  indicative  of 
mesoscale  flow  instabilities,  but  to  quantify  the 
role  of  flow  instabilities  required  further  analysis 
as  described  below. 

Metzger  and  Hurlburt  (2001)  describe  a  techni¬ 
que  to  separate  the  total  variability  field  of  model 
variables  into  two  components,  one  which  is  a 
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deterministic  response  to  the  forcing  and  the  other 
due  to  mcsoscale  flow  instabilities  (nondeterminis- 
tic).  The  method,  which  reveals  the  relative 
importance  of  the  two  anomaly-generation  me¬ 
chanisms,  is  summarized  as  follows.  First,  the 
long-term  mean  (1979  93)  SSH  at  each  gridpoint 
is  averaged  over  all  realizations.  This  mean  is  then 
subtracted  from  the  time  series  of  each  realization 
to  form  individual  time  series  of  anomaly  fields  (>/). 
At  each  mesh  point  the  temporally  corresponding 
records  of  all  realizations  are  averaged  (fj)  and 
deviations  (rf  =  >/  —  fj)  are  computed  about  the 
mean  comprised  of  all  realizations.  These  arc  then 
combined  and  averaged  as  follows: 

^  =  ?  +  ^  (4) 

where  the  bar  represents  the  average  overall 
realizations.  Averaging  over  the  same  time  period 
as  the  mean  formed  from  all  realizations,  as  was 
done  here,  ensures  that  the  temporal  average  of  the 
deterministic  anomalies  is  zero.  Now  the  right- 
hand  side  of  (4)  represents  the  deterministic 
variability  due  to  atmospheric  and  throughflow 
forcing  plus  the  variability  due  to  flow  instabilities, 
and  the  left-hand  side  is  the  total  variability. 
Subsequently,  each  component  can  be  divided  by 
the  total  variability  to  determine  the  fractional 
amount  of  variability  due  to  each,  allowing 
assessment  of  the  relative  impact  of  flow  instabil¬ 
ities  and  wind/th  rough  flow  forcing  on  model 
variability. 

As  explained  and  demonstrated  by  Metzger  and 
Hurlburt  (2001),  a  limited  number  of  realizations 
can  underestimate  the  nondeterministic  fraction 
of  variability  clue  to  flow  instabilities  and  over¬ 
estimate  the  deterministic  component.  Here,  we 
constructed  maps  of  the  percent  of  deterministic 
vs.  nondeterministic  variability  using  the  statistical 
technique  described  above  for  2  6  realizations. 
The  qualitative  difference  between  these  maps 
became  progressively  less  as  more  solutions  were 
included,  and  the  difference  between  including  five 
vs.  six  simulations  (not  shown)  was  quite  small. 

Fig.  1  1  shows  the  fractional  variability  for  SSH, 
surlacc-laycr  velocities,  and  abyssal-layer  eddy 
kinetic  energy  and  pressure.  The  high  values  (dark 
red)  indicate  regions  where  flow  instabilities  are 


the  primary  source  of  variability  and  the  dark 
blues  where  wind/through  flow  forcing  is  the 
dominant  source  of  variability.  In  Fig.  1  la  flow 
instabilities  are  the  dominant  source  of  SSH 
variability  over  much  of  the  subpolar  gyre  between 
37-38°N,  where  the  EKWC  separates  from  the 
coast  of  Korea,  and  following  the  general  pathway 
of  the  subpolar  front  along  the  north  rim  of  the 
Yamato  Rise.  The  nondeterministic  SSH  varia¬ 
bility  along  this  pathway  is  a  manifestation  of 
inertial  jets,  meandering,  and  eddy  activity  asso¬ 
ciated  with  the  subpolar  front.  The  result  that  the 
SSH  variability  in  the  region  where  the  EKWC 
separates  from  the  coast  of  Korea  is  primarily  due 
to  flow  instabilities  indicates  that  the  simulated 
variation  in  the  location  of  separation  (seasonal  or 
otherwise)  is  associated  mainly  with  instabilities  of 
the  jet  rather  than  seasonal  strengthening  (weak¬ 
ening)  of  the  subpolar  gyre  during  the  winter 
(summer)  in  response  to  wind  forcing  or  the 
seasonal  variation  in  the  straits  throughflow. 
However,  north  of  the  observed  separation  lati¬ 
tude  (>38°N)  the  SSH  variability  is  deterministic, 
suggesting  that  the  southward  extent  of  the 
subpolar  gyre,  especially  near  East  Korea  Bay,  is 
driven  by  the  wind  forcing  and  that  the  simula¬ 
tions  contain  little  or  no  eddy  shedding  from 
the  EKWC  after  it  separates  from  the  coast.  In 
Section  6,  we  found  that  northward  penetration 
of  the  post-separation  current  pathway  tends 
to  be  limited  by  an  eastward  abyssal  current 
along  38.5  38.8  N,  consistent  with  the  results  in 
Fig.  1  la-c. 

Other  regions  where  the  SSH  variability  is 
deterministic  include  the  deep  areas  over  the  Japan 
and  Yamato  Basins.  This  suggests  that  the  bottom 
topography  and  abyssal  flow  strongly  constrain 
the  surface  current  pathways  around  these  basins 
(Figs.  8  and  9).  The  ring  of  deterministic  SSH 
variability  around  the  edge  of  the  entire  basin  is 
due  to  wind-forced  coastally  trapped  internal 
Kelvin  waves  and  the  quasigeostrophic  constraint 
that  requires  constant  stream  function  along 
closed  boundary  segments  (Kelvin  waves  are  not 
quasigeostrophic).  The  dcterminicity  of  surface- 
layer  velocity  components  variability  (Figs,  lib 
and  c)  is  generally  similar  to  the  SSH  variability 
patterns,  except  in  the  region  of  the  UJleung  Basin. 
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Fig.  1 1  The  fraction  of  the  modeled  variability  from  four-layer  I  /32  NLOM  for  (A)  $SH.  (B)  //-component  of  surface  layer  velocity. 
(C)  /  -component  of  surface-layer  velocity.  (D)  abyssal-layer  eddy  kinetic  energy,  and  (E)  abyssal-layer  pressure  deviation.  The  blue 
areas  are  predominantly  deterministic  (caused  by  direct  wind  forcing)  and  the  red, orange  areas  are  predominantly  nondelerministic 
(caused  by  flow  instabilities). 


35N 


Immediately  surrounding  most  of  the  Ulleung 
Basin,  the  SSH  variability  is  highly  nondetermi- 
nistic.  Taken  with  the  pattern  of  nondeterminicity 
in  the  abyssal  layer,  this  pattern  is  consistent  with 
baroclinic  instability  surrounding  the  Ulleung 
Basin  in  the  surface  layer.  The  instability  transfers 
energy  to  the  abyssal  layer,  where  the  Plow  is 
inhibited  from  crossing  the  geostrophic  (///;) 
contours  of  the  topography  to  satisfy  conservation 


of  potential  vorticity.  Thus,  the  EKE  is  largely 
transferred  to  the  interior  of  the  Ulleung  Basin, 
and  the  resulting  KEM  (Fig.  9)  predominates  in 
following  the  jih  contours  of  the  bottom  topo¬ 
graphy.  The  relatively  uniform  pattern  of  abyssal 
pressure  nondeterminicity  covering  the  abyssal 
plain  of  the  Ulleung  Basin  and  nondeterminicity 
of  P4>  KE4  indicates  that  the  nondeterministic  P4 
variations  are  typically  the  scale  of  the  Ulleung 
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Basin  abyssal  plain.  Overall,  the  nondeterministic 
fraction  of  the  abyssal  variability  is  predominantly 
confined  to  areas  where  the  bottom  topography  is 
relatively  flat,  such  as  the  Japan  Basin  abyssal 
plain,  the  Yamato  Trough,  and  Ulleung  Basin. 

8.  Summary  and  conclusions 

Seven  different  monthly  wind-stress  climatolo¬ 
gies  (Hellcrman  and  Rosenstein,  COADS,  Na, 
ECMWF  10m,  ECMWF  lOOOmb,  NOGAPS  and 
NCEP)  were  used  to  force  a  suite  of  linear  and 
nonlinear  simulations  in  order  to  elucidate  the 
sensitivity  of  JES  simulations  to  the  choice  of 
wind- forcing  product,  including  the  impact  on 
the  model  dynamics.  Seasonally  varying  straits 
through  How  forcing  was  included  in  all  the 
simulations,  and  the  impact  of  throughflow 
forcing  in  the  absence  of  wands  was  also  investi¬ 
gated.  The  linear  solutions  provide  deterministic 
information  of  the  impact  of  the  w'ind  forcing 
without  the  complications  of  higher-order  dy¬ 
namics.  Those  higher-order  dynamics,  however, 
provide  the  ingredients  needed  for  the  realistic 
representation  of  many  of  the  circulation  features. 
Nonlinear  simulations  with  both  1/8C  resolution 
(w'cak  llow  instabilities)  and  1/32"  (strong  How 
instabilities)  w'erc  used. 

The  linear  simulations  include  only  the  lowest- 
ordcr  dynamics,  mainly  a  Munk  p1  1  viscous 
boundary  layer  with  a  Sverdrup  (1947)  interior. 
In  the  absence  of  wind  forcing,  the  TWC  hugs  the 
western  boundary  as  a  Munk  p1  3  viscous  bound¬ 
ary  layer  until  it  is  forced  to  separate  from  the 
coast  of  Korea  and  flow7  eastward  to  satisfy  the 
constraints  of  the  Tsugaru  and  Soya  outflow  ports. 
Thus,  contrary  to  what  is  observed,  there  is 
northeastward  llow  along  the  coast  of  Siberia 
and  extremely  weak  flow  along  the  coast  of 
Honshu.  When  w'inds  are  added,  in  all  cases  net 
positive  wind-stress  curl  in  the  northern  part  of 
the  domain  results  in  overall  cyclonic  flow  in  the 
subpolar  gyre,  including  southwestward  flow 
along  the  coast  of  Siberia.  In  most  cases  the 
separation  latitude  of  the  EKWC  is  still  deter¬ 
mined  by  the  latitude  of  Tsugaru  Strait  (too  far 
north),  although  some  of  this  flow  continues 


northward  on  the  eastern  side  of  the  basin  and 
exits  through  Soya.  The  same  simulations  also 
depict  southwestward  flow  along  part  or  all  of  the 
Honshu  coast,  contrary  to  observations.  The 
NCEP-forced  simulation  stands  alone  in  depicting 
northeastward  How  along  the  entire  Honshu  coast 
as  observed,  but  is  the  only  linear  simulation  that 
lacks  an  EKWC  and  overall  is  the  least  realistic  (a 
clear  outlier). 

The  nonlinear  simulations  include  the  addition 
of  bottom  topography,  multiple  internal  modes, 
diapycnal  mixing,  and  ventilation  of  layer  inter¬ 
faces.  At  l/8r  resolution,  these  additions  allow  the 
possibility  of  isopycnal  outcropping,  the  formation 
of  inertial  boundary  layers  and  jets,  and  weak 
barotropic/baroclinic  flow  instabilities.  At  this 
resolution,  mean  abyssal  currents  are  weak  and 
driven  predominantly  by  deep  vertical  mixing  in 
the  subpolar  gyre.  At  1/32  the  transfer  of  energy 
from  the  upper  ocean  to  the  abyssal  layer  via 
baroelinic  instability  is  much  stronger,  driving 
stronger  mean  abyssal  currents  that  have  a  more 
significant  influence  on  mean  surface  current 
pathways  than  at  1/8  ,  i.c.  much  stronger  upper 
ocean- topographical  coupling. 

In  general,  most  of  the  simulations  reproduce  the 
large-scale  circulation  features  of  the  JES  when 
compared  to  schematic  charts  of  the  circulation.  In 
particular,  all  depict  cyclonic  circulation  in  the 
subpolar  gyre  due  to  positive  wind-stress  curl, 
which  results  in  regional  upwelling  due  to  Ekman 
suction,  isopycnal  outcropping,  and  diapycnal 
mixing.  Likely  areas  of  intermediate-  and/or  deep¬ 
water  formation  occur  where  the  deepest  interface 
ventilates  (or  the  isopycnals  outcrop),  particularly 
over  a  long  time  period.  However,  substantial 
variability  in  the  strength,  size,  and  location  of 
the  subpolar  gyre  and  the  deep  interface  ventilation 
exists.  In  addition,  all  of  the  nonlinear  simulations 
(including  the  simulation  w'ith  straits  throughflow 
forcing  only)  simulate  the  northeastward  NB  of  the 
TWC  along  the  Honshu  coast.  This  is  due  to 
isopycnal  outcropping  in  the  simulations,  as 
explained  by  Hogan  and  Hurlburt  (2000). 

Simulation  of  the  EKWC  separation  latitude 
(37-38°N)  from  the  coast  of  Korea  has  proven  to 
be  one  of  the  more  challenging  problems  in  JES 
modeling.  Hogan  and  Hurlburt  (2000)  identified 
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two  distinct  mechanisms  that  could  yield  a  realistic 
separation  latitude  for  the  EKWC,  wind-stress 
cur]  and  upper  ocean  topographical  coupling,  the 
latter  only  occurring  in  the  1/32  simulations. 
However,  the  wind-stress  curl  mechanism  is 
ineffective  in  the  linear  simulations  and  only  one 
gives  a  separation  latitude  south  of  Tsugaru  Strait 
(—4 I  N).  That  is  the  NOGAPS-forced  simulation, 
which  gives  separation  at  ~39'  N.  In  this  study, 
we  found  that  when  part  of  the  flow  through 
Tsushima  Strait  is  diverted  into  the  NB  (as  in  all 
the  nonlinear  simulations),  the  wind-stress  curl 
mechanism  is  better  able  to  compete  with  the 
straits  throughflow  and  gives  a  realistic  EKWC 
separation  latitude  in  the  1/8  nonlinear  simula¬ 
tions  forced  by  Na,  ECMWF  10  m,  and  NOGAPS 
and  marginally  realistic  separation  with  ECMWF 
1000  mb  and  NCEP  forcing.  The  simulations 
forced  with  Hellerman-Rosenstein,  COADS.  and 
throughflow  only  still  show  EKWC  overshoot 
with  separation  at  ~40  N.  At  1/32  resolution  five 
of  the  six  simulations  show  realistic  EKWC 
separation,  the  Na,  ECMWF  10m,  and  NO¬ 
GAPS-forced  via  the  wind-stress  mechanism  as 
before,  the  Hellerman-Rosenstein  and  ECMWF 
1000  mb  via  the  upper  ocean  topographical  cou¬ 
pling  mechanism.  Hogan  and  Hurlburt  (2000) 
discuss  EKWC  separation  via  the  upper  ocean 
topographical  coupling  mechanism  for  the  Heller- 
man  Rosenstcin-forced  simulation  in  detail.  1/32 
throughflow  only  and  COADS-forced  simulations 
were  omitted  because  of  similarity  to  the  Heller- 
man  Rosenstein-forced  simulation  at  1/8  resolu¬ 
tion.  These  would  probably  also  give  a  realistic 
EKWC  separation  latitude  at  1/32  resolution  via 
the  upper  ocean- topographical  coupling  mechan¬ 
ism.  At  1/32  resolution  only  the  NCEP-forced 
simulation  shows  EKWC  overshoot  (to  ^43  N). 
This  is  due  to  strong  wind-stress  curl  support  for  a 
northward  western  boundary  current  between 
39rN  and  43  N  as  seen  in  the  linear  simulation. 

At  l/32e  all  of  the  simulations  except  the  NCEP- 
forced  show  a  strong  influence  of  the  Yamato 
Rise  on  the  path  of  the  subpolar  front  due  to 
steering  by  abyssal  currents.  The  1/32  NCEP- 
forced  simulation  docs  not  show  this  because  the 
subpolar  front  approaches  the  Yamato  Rise  from 
the  northwest  in  that  case.  The  path  of  the 


subpolar  front  in  the  1/8  NOGAPS  and  ECMWF 
lOm-forced  simulations  is  also  influenced  by 
the  Yamato  Rise.  Of  the  1/8  simulations,  the 
NOGAPS-forced  simulation  has  the  strongest 
abyssal  currents.  These  are  driven  by  the  strongest 
isopycnal  outcropping  and  vertical  mixing  in  the 
subpolar  gyre.  In  this  case,  some  of  the  vorticity 
produced  by  the  wind-stress  curl  is  transferred 
directly  into  the  abyssal  layer.  This  is  the  reason  it 
is  quite  similar  to  the  corresponding  l/32c  simula¬ 
tion  despite  differences  in  the  abyssal  forcing 
mechanisms  (isopycnal  outcropping  at  1/8  and 
strong  baroclinic  instability  at  1/32  ).  In  the  1/8' 
ECMWF  lOm-forced  simulation,  weaker  abyssal 
currents  around  the  Yamato  Rise  are  aided  by  an 
overlying  region  of  strong  negative  wind-stress 
curl  not  seen  in  the  other  wind  sets. 

Synoptic  wind  forcing  was  used  to  investigate 
determinicity  vs.  nondeterminicity  of  interannual 
variability  in  1/32  simulations  forced  with  the 
1979-93  six-hourly  ECMWF  10-m  reanalysis  wind- 
stresses  and  seasonally  varying  flows  through  the 
straits.  Mean  SSH  from  six  simulations  that 
differed  only  in  the  initial  state  show  significantly 
different  mean  current  pathways  in  1-year  means 
over  the  same  year.  The  differences  are  primarily 
due  to  flow  instabilities,  which  are  well  resolved  at 
1/32  resolution.  The  fraction  of  variability  due  to 
flow  instabilities  (the  nondeterministic  component) 
and  wind/throughflow  forcing  (the  deterministic 
component)  was  calculated,  and  indicates  that  flow 
instabilities  are  responsible  for  most  of  the  SSH 
variability  in  the  region  where  the  EKWC  separates 
from  the  coast  of  Korea  (^ — 37  -38  N),  extending 
along  the  pathway  of  the  subpolar  front  to  the 
Tsugaru  Strait.  These  instabilities  are  associated 
with  the  inertial  character  of  the  EKWC  recircula¬ 
tion  gyre  (the  Ulleung  eddy)  and  meandering  and 
frontal  eddy  activity  associated  with  the  subpolar 
front.  SSH  variability  is  largely  deterministic  in  the 
coastal  Kelvin  waveguide  and  over  some  deep 
basins  and  troughs.  Abyssal  variability  of  pressure 
and  eddy  kinetic  energy  was  nondeterministic  over 
areas  with  relatively  flat  topography,  indicative  of 
the  constraint  that  deep  flows  tend  to  follow 
geostrophic  (///?)  contours  when  |V(///?)|  is  large. 

Collectively,  the  simulations  forced  with  the 
different  wind  sets  show  a  much  wider  range  of 


1488 


P.J.  Hogan,  H.E.  Huvlburi  /  Deep-Sea  Research  II  52  (2005)  1464-1489 


features  and  dynamics  than  obtained  using  a  single 
wind  set.  Thus,  they  give  a  broader  perspective 
on  the  JES  circulation  and  the  range  of  features 
and  dynamics  it  may  exhibit  than  found  in  the 
simulations  driven  by  any  single  atmospheric 
forcing  product.  Based  on  the  mean  circulation, 
overall  the  most  realistic  simulations  (all  nonlinear) 
are  the  1/32°  ECMWF  1000  mb  forced,  the  1/8° 
and  1/32°  Na  forced,  and  the  1/32°  NOGAPS 
forced.  The  least  realistic  nonlinear  simulations  are 
the  1/8°  and  1/32°  NCEP  forced.  This  is  based  on 
general  agreement  with  the  schematics  and  other 
known  features  of  the  circulation.  Although  wind 
forcing  strongly  affects  the  circulation  features 
and  current  pathways  in  the  JES,  other  mechan¬ 
isms,  both  physical  and  dynamical,  also  play  an 
important  role.  The  upper  ocean-topographical 
coupling  described  in  Hogan  and  Hurlburt  (2000) 
is  an  example  of  a  dynamical  mechanism.  Other 
physical  mechanisms  such  as  buoyancy  forcing, 
accurate  mixed-layer  representation,  and  realistic 
straits  throughflow,  and  the  nonlinear  combination 
of  all,  undoubtedly  play  an  important  role  as  well. 
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